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Abstract 

Usher syndrome is an autosomal recessive disorder manifesting hearing loss, retinitis pigmentosa and vestibular 
dysfunction, and having three clinical subtypes. Usher syndrome type 1 is the most severe subtype due to its profound 
hearing loss, lack of vestibular responses, and retinitis pigmentosa that appears in prepuberty. Six of the corresponding 
genes have been identified, making early diagnosis through DNA testing possible, with many immediate and several long- 
term advantages for patients and their families. However, the conventional genetic techniques, such as direct sequence 
analysis, are both time-consuming and expensive. Targeted exon sequencing of selected genes using the massively parallel 
DNA sequencing technology will potentially enable us to systematically tackle previously intractable monogenic disorders 
and improve molecular diagnosis. Using this technique combined with direct sequence analysis, we screened 1 7 unrelated 
Usher syndrome type 1 patients and detected probable pathogenic variants in the 16 of them (94.1%) who carried at least 
one mutation. Seven patients had the MY07A mutation (41.2%), which is the most common type in Japanese. Most of the 
mutations were detected by only the massively parallel DNA sequencing. We report here four patients, who had probable 
pathogenic mutations in two different Usher syndrome type 1 genes, and one case of MY07A/PCDH15 digenic 
inheritance. This is the first report of Usher syndrome mutation analysis using massively parallel DNA sequencing and the 
frequency of Usher syndrome type 1 genes in Japanese. Mutation screening using this technique has the power to quickly 
identify mutations of many causative genes while maintaining cost-benefit performance. In addition, the simultaneous 
mutation analysis of large numbers of genes is useful for detecting mutations in different genes that are possibly disease 
modifiers or of digenic inheritance. 
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Introduction 

Usher syndrome (USH) is an autosomal recessive disorder 
characterized by hearing loss (HL), retinitis pigmentosa (RP) and 
vestibular dysfunction. Three clinical subtypes can be distin- 
guished. USH type 1 (USHl) is the most severe among them 
because of profound HL, absent vestibular responses, and 
prepubertal onset RP. USH type 2 (USH2) is characterized by 
congenital moderate to severe HL, with a high-frequency sloping 
configuration. The vestibular function is normal and onset of RP is 



in the first or second decade. The onset of the visual symptoms 
such as night blindness in USH usually occurs several years later 
than in USHl. USH type 3 (USH3) is characterized by variable 
onset of progressive HL, variable onset of RP, and variable 
impairment of vestibular function (normal to absent) [1,2]. 

To date, nine genetic loci for USHK^USHIB-H, J, and A) have 
been mapped to chromosomes llql3.5, llpl5.1, 10q22.1, 21q21, 
10q21-q22, 17q24-q25, 15q22-q23 {USHIH and J), and 
lOpl 1.21-q21.I [2,3,4]. Six of the corresponding genes have 
been identified: the actin-based motor protein myosin Vila 



PLOS ONE I www.plosone.org 



1 



March 2014 | Volume 9 | Issue 3 | e90688 



MPS for Diagnosis of Usher Syndrome Type 1 



(MT07A, USHIB) [5]; two cadherin-related proteins, cadherin 23 
{CDH23, USHID) [6] and protocadherin 15 {PCDH15, USHIF) 
[7]; and two scaffold proteins, harmonin [USHIC) [8] and sans 
{USHiq [9]; the Ca^^- and integrin-binding protein {CIB2, 
USHIJ) [4]. In Caucasian USHl patients, previous studies 
showed that mutations in Mr07A, USHIC, CDH23, PCDH15, 
and USHIG, were found in 39-55%, 7-14%, 7-35%, 7-1 1 %, and 
0-7%, respectively (the frequency of CIB2 is still unknown) 
[10, 11, 12]. In Japanese, Nakanishi et al. showed that MY07A and 
CDII23 mutations are present in USHl patients [13], however, 
the frcquc'ncy is not yet known. In addition, mutations in three 
corresponding genes (usherin USH2A [14], G protein-coupled 
receptor 98; GPR98 [15], and deafness, autosomal recessive 31; 
DFNB31 [16]) have been reported so far in USH2, and USH3 is 
caused by mutations in the clarin 1 (CLRNT) [17] gene. 

Comprehensive molecular diagnosis of USH has been ham- 
pered both by genetic heterogeneity and the large number of 
exons for most of the USH genes. The six USHl genes collectively 
contain 180 coding exons [4,9,10] the three USH2 genes comprise 
175 coding exons [15,16,18], and the USH3 gene has five coding 
exons [17]. In addition some of these genes are alternatively 
spliced ([4,7,8,16,17] and NCBI database: http://www.ncbi.nlm. 
nih.gov/nuccore/). Thus far, large-scale mutation screening has 
been performed using direct sequence analysis, but that is both 
time-consuming and expensive. We thought that targeted exon 
sequencing of selected genes using the Massively Parallel DNA 
Sequencing (MPS) technology would enable us to systematically 
tackle previously intractable monogenic disorders and improve 
molecular diagnosis. 

Therefore, in this study, we have conducted genetic analysis 
using MPS-based genetic screening to find mutations in nine 
causative USH genes (except CIB2) in Japanese USHl patients. 

Results 

Mutation analysis of the nine USH genes in 1 7 unrelated USH 1 
patients revealed 19 different probable pathogenic variants, of 
which 14 were novel (Table 1). 

AU mutations were detected in only one patient each and 
sixteen of the 1 7 patients (94. 1 %) carried at least one mutation, 
while one patient had no mutations. Thirteen of the 1 6 mutation 
carriers each had two patliog(;nic mutations (Table 2). 

Nonsense, frame shift, and splice site mutations are all classified 
as pathogenic, whereas missense mutations are presumed to be 
probable pathogenic variants based on results of prediction 
software for evaluation of the pathogenicity of missense variants 
(Table 1). 

Of the 19 probable pathogenic mutations that we found, 17 
were detected by MPS. The remaining two (p.Lys542GlnfsX5 in 
Mr07A and c.5821-2A>G in CDH23) were sequenced by direct 
sequence analysis. 

Of our 17 USH patients, seven had Mr07A mutations (41.2%), 
tiiree had CDH23 mutations (17.6%), and two had PCDH15 
mutations (11.8%). We did not find any probable pathogenic 
mutations in USHIC, USHIG, and USH2/3 genes. 

Four USHl patients (Cases #3, 5, 8, 15) had probable 
pathogenic mutations in two different USH genes, with one being 
a biallelic mutation (Table 3). The other heterozygous/homcjzy- 
gous mutations were missense variants. Three of these patients 
(Cases #3, 5, 8) presented with earlier RP onset (night blindness) 
than in the other patients with two pathogenic mutations (Cases 
#1, 6, 7, 9, 11, 16) (/» = 0.007) (Fig. 1). 

One patient (Case #4) had heterozygote mutations in two 
USHl genes (p.Ala771Ser in Mr07A and C.158-1G>A in 



PCDHIS). His parents and one brother were found to also be 
carriers for these mutations. Another brother had no variants 
(Fig. 2). 

Discussion 

For USHl, early diagnosis has many immediate and several 
long-t(;rm advantages for patients and their families [1]. However, 
diagnosis in childhood, based on a clinical phenotype, has been 
difficult because patients appear to have only non-syndromic HL 
in childhood and RP develops in later years. Although early 
diagnosis is now possible through DNA testing, performing large- 
scale mutation screening for USH genes in all non-syndromic HL 
children has been both time-consuming and expensive. Therefore, 
the availability of MPS, which facilitates comprehensive large- 
scale mutation screening [19] is a very welcome advance. 

MPS technology' enabled us to detect pathogenic mutations in 
USHl patients efficiently, identifying one or two pathogenic/likely 
pathogenic mutations in 16 of 17 (94.1%) cases. This was 
comparable to previous direct sequence analysis results such as 
Bonnet et al. who detected one or two mutations in 24 out of 27 
(89'1'i)) USHl patients [11] and Le Quesne Stabej et al. who 
detected one or two mutations in 41 out of 47 (87.2%) USHl 
patients [12]. 

In addition, MPS assists in the analysis of disease modifiers and 
digenic inheritance because it simultaneously investigates many 
causative genes for a specific disease, such as in our case, USH. 
Previous reports have described several USH cases with patho- 
genic mutations in two or three different USH genes [11,1 2,20] . In 
our study, four patients had two pathogenic mutations in one gene 
and missense variants in a different gene (Table 3). We considered 
the latter to possibly be a disease modifier. For example, 
L'S//iC:p.Tyr813Asp, which occurred in 0/384 control chromo- 
somes and was predicted to be "probably damaging" by the 
Poly-phen program, was found with a homozygous CDH23 
nonsense mutation (p.Arg2107X) (Case #15). As for what the 
variant "modifies", we speculate that for USHl patients with a 
disease modifier, RP symptoms such as night bUndness show an 
earlier onset. However, we think that profound HL and the 
absence of vestibular function in USH 1 patients are not affected 
by modifiers as they are congenital and therefore not progressive. 

Ebermann et al. described a USH2 patient with "digenic 
inheritance." a heterozygous truncating mutation in GPR98, and a 
truncating heterozygous mutation in PDZ domain-containing 7 
{PDZDT}, which is reported to be a cause of USH [20]. Our USHl 
patient (Case #4) had segregated MrO 7^:p.Ala7 7 ISer and 
PCD//75:c. 158-lOA. Molecular analyses in mouse models have 
shown many interactions among the USHl proteins [2]. In 
particular, MT07A directiy binds to PCDHIS and both proteins 
are expressed in an overlapping pattern in hair bundles in a mouse 
model [21]. PCDffl5:c.l58-lG>A, predicted to alter the splice 
donor site of intron 3, has been classified as pathogenic. 
MTOZ4:p.Ala771Ser is a non-truncating mutation, but was 
previously reported as disease-causing [13]. So, we consider the 
patient to be the first reported case ot MY07A/ PCDHIS digenic 
inheritance. 

However, we should be aware of two limitations of MPS 
technology. First, the target region of MPS cannot cover all coding 
exons of USH genes. Actually, the coverage of the target exons 
was 97.0% in our study. So, it is impossible to detect a mutation in 
a region which is not covered using this system (Case #9: c.5821- 
2A>G). Secondarily, the MPS system used in this study, is not 
effective for detecting homo-polymer regions, for example poly C 
stretch [22] (Case #8: p.Lys542GlnfsX5). In addition, concerning 
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Table 2. Details of phenotype and genotype of 17 USHl patients 







Sample No. 


Age 


Sex 


Allelel 


Allele2 


Hereditary 
form 


Onset of 

night 

blindness 


Cataract 


Cochlear 
Hearing Aid Implant 


MY07A 


1 


37 


M 


p.Gln493X 


p.Trpl495CysfsX55 


sporadic 


13 


no 


unilateral 


unilateral 


2 


41 


W 


p.l2069fsX6 


p.l2069fsX6 


AR 


unknown 


both eyes 


bilateral 


no 


5 


54 


M 


p.Val692Met 


p.Val692Met 


AR 


5 


both eyes 


no 


no 


6 


54 


W 


p.Arg570X 


p.Arg570X 


sporadic 


6 


no 


no 


no 


8 


14 


M 


p.Lys542GlnfsX5 


p.Lys542GlnfsX5 


sporadic 


6 


no 


unilateral 


unilateral 


n 


54 


M 


p.Asp2010Asn 


p.Trp2107X 


sporadic 


13 


no 


no 


no 


17 


56 


W 


p.Cys705X 


p.Cys705X 


sporadic 


unknown 


no 


no 


no 


CDH23 


7 


12 


W 


p.Argl189ArglfsX5 


p.Argl189ArglfsX5 


sporadic 


12 


both eyes 


no 


bilateral 


9 


9 


M 


p.Ser1927Cysfs16 


c.5821-2A>G 


sporadic 


8 


no 


unilateral 


unilateral 


15 


16 


W 


p.Arg2107X 


p.Arg2107X 


sporadic 


unknown 


no 


no 


no 


PCDH15 


3 


47 


W 


p.Glul113X 


p.Glu1113X 


sporadic 


5 


both eyes 


no 


no 


16 


28 


W 


p.Arg991X 


p.Arg991X 


AR 


10 


no 


no 


no 


10 


62 


M 


p.Arg962Cys 


unknown 


sporadic 


9 


both eyes 


no 


no 


12 


52 


M 


p.Arg336X 


unknown 


sporadic 


3 


no 


no 


no 


13 


51 


M 


p.Val1242Met 


unknown 


sporadic 


10 


no 


no 


no 




4 


21 


M 


p.Ala771Ser»' 


C.I 58-1 G>A«^ 


sporadic 


10 


no 


unilateral 


unilateral 


unknown 


14 


64 


W 


unknown 


unknown 


sporadic 


15 


both eyes 


unilateral 


no 



*AII subjects have congenital deafness and RP. 
doi:l 0.1 371 /journal.pone.0090688.t002 



pathogenecity of mutations identified, functional analysis will be 
necessary to draw the final conclusion in the future. 

In UK and US Caucasian USHl patients, USHIB [MY07A) 
has been reported as the most common USHl genetic subtype 
[11,12], while USHIF {PCDH15) has been reported as the most 
common USHl genetic subtype in North American Ashkenazi 
Jews [23]. In Japanese, our study revealed that the most common 
type was MY07A (41.7%), which was similar to the frequency in 
the above Caucasian patients (46.8~55%) [11,12]. However, the 
small number of USHl patients in our study might have biased the 
frequency and further large cohort study will be needed in the 
future. 

In addition, most of our detected mutations were novel. We 
have previously reported genes responsible for deafness in 
Japanese patients and observed differences in mutation spectrum 
between Japanese (who are probably representative of other Asian 
populations) and populations with European ancestry [24] . 

In conclusion, our study was the first report of USH mutation 
analysis using MPS and the frequency of USHl genes in Japanese. 
Mutation screening using MPS has the potential power to quickly 
identify mutations of many causative genes such as USH while 
maintaining cost-benefit performance. In addition, the simulta- 
neous mutation analysis of large numbers of genes was useful for 
detecting mutations in difierent genes that are possibly disease 
modifiers or of digenic inheritance. 



Materials and Methods 

Subjects 

We screened 17 Japanese USHl patients (aged 9 to 64 years): 
three from autosomal recessive families (non-aflfected parents and 
two or more affected siblings), and 14 from sporadic families. 
There were 9 males and 8 females. None of the subjects had any 
other noteworthy symptoms. AU subjects or next of kin on the 
behalf of the minors/children gave prior written informed consent 
for participation in the project, and the Ethical Committee of 
Shinshu University approved the study and the consent procedure. 

Amplicon Library Preparation 

An Amplicon library of the target exons was prepared with an 
Ion AmpliSeq Custom Panel (Applied Biosystems, Life Technol- 
ogies, Carlsbad, CA) designed with Ion AmpliSeq Designer 
(https://www.ampliseq.com/browse.action) for nine USH genes 
by using Ion AmpliSeq Library Kit 2.0 (Applied Biosystems, Life 
Technologies) and Ion Xpress Barcode Adapter 1-16 Kit (Applied 
Biosystems, Life Technologies) according to the manufacturers' 
procedures. 

In brief, DNA concentration was measured with Quant-iT 
dsDNA HS Assay (Invitrogen, Life Technologies) and Qubit 
Fluorometer (Invitrogen, Life Technologies) and DNA quality was 
confirmed by agarose gel electrophoresis. 10 ng of each genomic 
DNA sample was amplified, using Ion AmpliSeq HiFi Master Mix 
(Applied Biosystems, Life Technologies) and AmpliSeq Custom 
primer pools, for 2 min at 99°C, followed by 15 two-step cycles of 
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Figure 1. The number of mutations and the age of RP onset in 
Usher syndrome type 1 patients. The age of RP onset is earlier in 
the patients with more than two pathogenic mutations. RP: retinitis 
pigmentosa. 

doi:10.1371/journal.pone.0090688.g001 
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99°C for 15 sec and 60°C for 4 min, ending with a holding period 
at 10°C in a PGR thermal cycler (Takara, Shiga, Japan). After the 
Multiplex PGR amplification, amplified DNA .samples were 
digested with FuPa enzyme at 50°C for 10 min and 55°G for 
10 min and the enzyme was successively inactivated for 60°C for 
20 min incubation. After digestion, diluted barcode adapter mix 
including Ion Xpress Barcode Adapter and Ion PI adaptor were 
ligated to the end of the digested amplicons with ligase in the kit 
for 30 min at 22°C and the ligase was successively inactivated at 
60°C for 20 min incubation. Adaptor ligated amplicon libraries 
were purified with the Agencourt AMPure XP system (Beckman 
Coulter Genomics, Danvers, MA). The amplicon libraries were 
quantified by using Ion Library Quantitation Kit (Applied 
Biosystems, Life Technologies) and the StepOne plus realtime 
PGR system (Applied Biosystems, Life Technologies) according to 
the manufacturers' procedures. After quantification, each ampli- 
con library was diluted to 20 pM and the same amount of the 1 2 
libraries for 1 2 patients were pooled for one sequence reaction. 

Emulsion PCR and Sequencing 

The emulsion PGR was carried out with the Ion OneTouch 
System and Ion OneTouch 200 Template Kit v2 (Life Technol- 
ogies) according to the manufacturer's procedure (Publication Part 
Number 4478371 Rev. B Revision Date 13 June 2012). After the 
emulsion PCR, template-positive Ion Sphere Particles were 
enriched with the Dynabeads MyOne Streptavidin Gl Beads 
(Life Technologies) and washed with Ion OneTouch Wash 
Solution in the kit. This process were performed using an Ion 
OneTouch ES system (Life Technologies). 

After the Ion Sphere Particle preparation, MPS was performed 
with an Ion Torrent Personal Genome Machine (PGM) system 
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Figure 2. Pedigree and sequence chromatograms of the patient witli tKie p.Ala771Ser in MY07A and C.158-10A in PCDH15 
mutations. (A) The pedigree and sequence results of the proband and family. (B) Sequence chromatograms from wild-type and mutations. The 
proband, his mothor and one brother carried a heterozygous 231 1G>T transition in exon 20, which results in an alanine to a serine (Ala771Ser) in 
MY07A. Another variation, 158-10A in intron 3 of PCDHIS, was derived from the proband and his father. Another brother had no variants. 
doi:1 0.1 371 /journal.pone.0090688.g002 



using the Ion PGM 200 Sequencing Kit and Ion 318 Chip (Life 
Technologies) according to the manufacturer's procedures. 

Base Call and Data Analysis 

The sequence data were processed witli standard Ion Torrent 
Suite Software and Torrent Server successively mapped to human 
genome sequence (buUd GRCh37/hgl9) with Torrent Mapping 
Alignment Program optimized to Ion Torrent data. The average 
of 562.33 Mb sequences with about 4,300,000 reads was obtained 
by one Ion 318 chip. The 98.0% sequences were mapped to the 
human genome and 94% of them were on the target region. 
Average coverage of depth in the target region was 314.2 and 
93.8% of them were over 20 coverage. 

After the sequence mapping, the DNA variant regions were 
piled up with Torrent Variant Caller plug-in software. Selected 
variant candidates were filtered with the average base QV 
(minimum average base quality 25), variant frequency (40-60% 
for heterozygous mutations and 80-100% for homozygous 
mutations) and coverage of depth (minimum coverage of depth 
10). After the filtrations, variant effects were analyzed with the 
wANNOVAR web site [25,26] (http://wannovar.usc.edu) includ- 
ing the functional prediction software for missense variants: 
Sorting Intolerant from Tolerant (SIFT; http://sift.jcvi.org/), 
and Polymorphism Phenotyping (PolyPhen2; http://genetics.bwh. 
harvard.edu/pph2/). The sequencing data was available in the 
DNA databank of Japan (Accession number: DRA001273). 



Algorithm 

Missense, nonsense, and splicing variants were selected among 
the identified variants. Variants were further selected as less than 
1% of 1) the 1000 genome database (http://www. lOOOgenomes. 
org/), 2) the 5400 exome variants (http://evs.gs.washington.edu/ 
EVS/), and 3) the in-house control. Candidate mutations were 
confirmed by Sanger sequencing and the responsible mutations 
were identified by segregation analysis using samples from family 
members of the patients. In addition, the cases with heterozygous 
or no causative mutation were fully sequenced by Sanger 
sequencing for USHl genes in order to verify the MPS results. 

Direct Sequence Analysis 

Primers were designed with the Primer 3 plus web server 
(http:/ /www.bioinformatics.nl/cgi-bin/primer3plus/primer3plus. 
cgi). Each genomic DNA sample (40 ng) was amplified using 
Ampli Taq Gold (Life Technologies) for 5 min at 94°C, followed 
by 30 three-step cycles of 94°C for 30 sec, 60°C for 30 sec, and 
72°C for 30 sec, with a final extension at 72°C for 5 min, ending 
with a holding period at 4°C in a PGR thermal cycler (Takara, 
Shiga, Japan). The PGR products were treated with ExoSAP-IT 
(GE Healthcare Bio, Buckinghamshire, UK) and by incubation at 
37°C for 60 min, and inactivation at 80°C for 15 min. After the 
products were purified, we performed standard cycle sequencing 
reaction with ABI Big Dye terminators in an ABI 3130x1 
sequencer (Life Technologies). 
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Accession numbers 

Mr07A, [NM_000260.3]; USHIC, [NM_ 153676.3]; CDH23, 
[NM_ 022124.5]; PCDH15, [NM_ 033056.3]; USHIG, [NM. 
173477.2]; USH2A, [NM_206933.2]; GPR98, [NM_ 032119.3]; 
DFNB31, [NM_ 015404.3]; CLRNl, [NM_ 174878.2]; PDZD7, 
[NM_ 001195263.1]. 
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